Hindered mobility of a particle near a soft interface 



Thomas BickeQ 
CPMOH, Umversite Bordeaux 1 & CNRS (UMR 5798) 
351 cours de la Liberation, 33405 Talence, France 



O 
O 

(N 
-i— > 
O 

O 

(N 



o 

CO 



o 
o 



O 

o 

-I— > 
c3 



i 

-a 
c 

o 
o 



X 



The translational motion of a solid sphere near a deformable fluid interface is studied in the 
low Reynolds number regime. In this problem, the fluid flow driven by the sphere is dynamically 
coupled the instantaneous conformation of the interface. Using a two-dimensional Fourier transform 
technique, we are able to account for the multiple backflows scattered from the interface. The 
mobility tensor is then obtained from the matrix elements of the relevant Green function. This 
analysis allows us to express the explicit position and frequency dependence of the mobility. We 
recover in the steady limit the result for a sphere near a perfectly flat interface. At intermediate time 
scales, the mobility exhibits an imaginary part, which is a signature of the elastic response of the 
interface. In the short time limit, we find the intriguing feature that the perpendicular mobility may, 
under some circumstances, become lower than the bulk value. All those results can be explained 
from the definition of the relaxation time of the soft interface. 



PACS numbers: 82.70.Dd, 68.05.-n, 47.15. G- 



I. INTRODUCTION 

The motion of a particle in the vicinity of a bounding 
surface is a long standing problem in colloidal science . 
When a colloidal sphere suspended in a quiescent fluid 
approaches a wall, the drag force acting on it increases 
with respect to the drag force when far from the wall. 
This property is attributed to hydrodynamic interactions 
that develop because of the boundary conditions imposed 
by the wall on the fluid flow. In addition, the motion 
of the particle becomes anisotropic since the mobility is 
higher in the direction parallel to the wall than in the 
perpendicular direction. 

Although the first investigations on the influence of a 
bounding wall date back to the early work of Lorentz Q , 
this field has known a certain revival during the past 
two decades. The main reason for this is certainly the 
achievement of technical progress, in particular in the 
field of single-molecule techniques, that allows nowa- 
days to measure the position-dependent mobility of in- 
dividual micron-size particles with a great accuracy. 
Among the most efficient tools, one can quote evanes- 
cent waves techniques 0, 0, single particle tracking 
by video-microscopy p| , particle handling with optical 
tweezers 
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noise analysis [9|, or fluores- 
cent correlation spectroscopy [ljj. Those various meth- 
ods share the common feature of probing the random 
motion of Brownian objects near one or two solid wall, 
and the the mobility coefficients deduced from the exper- 
imental data agree remarkably well with the theoretical 
predictions Q. 

The renewal of interest for this question is also due 
to the development of microfluidics Indeed, by re- 
ducing the size of the systems, the influence of surface 
effects are inevitably enhanced with respect to the bulk 



'Electronic address: th.bickel@cpmoh.u-bordeauxl.fr 



properties. Consequently, most of the physical phenom- 
ena take place near the boundaries. A fundamental un- 
derstanding on how surface properties might affect the 
overall flow field has therefore become crucial in order to 
propose new solutions that would take advantage of this 
predominance. Lastly, it has been suggested recently to 
use colloidal particles as local probes of the flow proper- 
ties near surfaces. This idea has been introduced in the 
context of the no-slip boundary condition [l2( , where the 
motion of the particles is expected to contain a signa- 
ture of the slip length 0, . More generally, one can 
think of a Brownian particle as a probe of the viscoelastic 
properties of the bounding surface. 

From a theoretical viewpoint, the motion of a solid par- 
ticle in the presence of a nearby, plane interface has been 
extensively studied in the past. During the last few years, 
calculations of mobility coefficients have been extended 
to particles near surfactant-covered interfaces [lif. in a 
liquid film between two fluids 15], or in a Poisseuile flow 
between planar walls 1161 . The effect of fluid inertia has 
also been accounted for Il7l. a s well as the possibility of 
liquid slippage at the wall Here, we re-examine this 
question for a particle near a fluid- fluid interface. Re- 
sults are available for the drag force acting on a sphere 
of radius a moving at a distance zq of a perfectly flat in- 
terface, up to second order in the ratio cl/zq [T||. While 
this problem is of some intrinsic interest, and is a logical 
starting point in the limit of very high surface tension, 
it is obvious that a real interface will generally deform 
owing to the motion of the particle. For finite surface 
tension, the motion of the particle is expected to be dy- 
namically coupled to the conformations of the interface. 
Indeed, the fluid flow caused by the displacement of the 
particle exerts stresses that deform the interface. Relax- 
ing back to its equilibrium position, the interface creates 
a backflow that in turn perturbs the motion of the par- 
ticle, and so on. The delay in the response of the soft 
surface to hydrodynamic stresses is therefore expected to 
induce memory effects in the motion of the particle |19| . 
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In general, the problem of the motion near a soft 
surface is highly non-linear due to the fact that the 
shape of the interface is unknown. Althought it cannot 
be solved exactly, iterative solutions have been derived 
when the deformation of the interface is asymptotically 
small [2£j, El- The idea is to first solve the motion of 
a spherical bead near a flat surface. As the resulting 
velocity produces an imbalance of normal stress at the 
interface, it is then possible to determine a first nonzero 
approximation for the deformation po| . This strategy is 
however limited as it only describes the first "image" cor- 
rection to hydrodynamic interactions. Also, it assumes a 
quasi-steady deformation profile and does not allow for 
a possible delay inherent in the response of an elastic 
interface. 

In this article, we present an analytical method that 
rigorously accounts for the infinite series of hydrody- 
namic reflections on the soft interface. This scheme is 
achieved within the only assumption that interface de- 
formations remain moderate. The rest of the paper is 
organized as follows. In Section^ we specify the system 
and introduce the general set of equations that govern the 
problem. We reformulate in Section ITTT1 the small defor- 
mation problem in terms of equivalent boundary condi- 
tions at the undisplaced interface. Results for the Oseen 
tensor and the mobility coefficients are then discussed in 
Section IIVI In particular, we find that the frequency- 
dependent mobility switches between two regimes over 
a time scale corresponding to the relaxation time of the 
interface. Finally, we come back to the relationship with 
experiments and draw some concluding remarks in Sec- 
tion El 



II. FORMULATION OF THE PROBLEM 
A. Linear hydrodynamics 

We consider a spherical particle of radius a moving 
near a fluid interface in the low Reynolds number regime. 
The interface separates two viscous, incompressible and 
immiscible fluids. Its average position is chosen to coin- 
cide with the x - y plane, with the z-coordinate directed 
perpendicular to the it. The two fluids are labelled with 
indices 1 and 2, fluid 1 lying above fluid 2. Further- 
more, we denote 771 and 772 the shear viscosities, p\ and 
pi the mass densities, and Ap = P2 — P\ > the mass 
density difference. In order to get the mobility tensor of 
the particle, we shall first evaluate the appropriate Green 
function — called the Oseen tensor in this context — and 
investigate the effect of a time-dependent point force F(t ) 
acting at position ro = (xo,yo,zo) on the flow field |22| . 
Without loss of generality, we can assume that the sphere 
is fully immersed in fluid 1. For small-amplitude and 
low-frequency motion, the flow velocity v(r, t) and the 



pressure p(r, t) are governed by the Stokes equations 

r/ Q V 2 v - Vp + F5{r - r ) = , (1) 
V.v = 0, (2) 

with a = 1 or 2, depending on whether the point r is 
located above or below the interface. In Eq. (Q, S stands 
for the Dirac delta-distribution. The two fluids are as- 
sumed to be quiescent except for the disturbance flow 
caused by the motion of the sphere. 

B. Physics of interfaces 

The Stokes equations have to be solved together with 
the usual boundary conditions at the interface, namely 
the velocity and the tangential constraints must be con- 
tinuous. The normal- normal component of the stress 
tensor presents a discontinuity which is balanced by the 
restoring force exerted by the deformed interface on the 
fluid. This question is quite involved since, in principle, 
the tangential and normal directions depend on the local 
and instantaneous conformation of the interface. How- 
ever, an approximate solution can be found for moderate 
deformations. In this case, the position of the almost 
flat interface can be described by a single-valued func- 
tion h(p,i), with p = {x,y). For our purpose, it is more 
convenient to use the two-dimensional Fourier represen- 
tation 

h(q,t) = J d 2 pexp{-iq.p]h(p,t) , (3) 

with q = (q x , q y ). The elastic properties of the interface 
are then described by the Hamiltonian |2^| 

U=\ |d 2 q( g 2 + Z- 2 )|Mq,*)f, (4) 

where 7 is the surface tension and l c — y/j/JgAp) the 
capillary length, g being the gravitational acceleration. 
The capillary length scale typically lies in the millimeter 
range for 7 w 100 mN/m, but can be as low as a few 
microns for ultra-soft interfaces with 7 w 0.1 pN/m [24j. 
We then proceed in the same manner as for linearized 
theory of capillary waves and express all the boundary 
conditions at the undisplaced interface = 0. This hy- 
pothesis of smooth deformation is valid up to linear order 
in the deformation field h, so that our approach is fully 
consistent with the harmonic description of the interface 
energy Eq. Q. 

C. Method of solution 

In spite of these classical simplifications, the coupling 
between the motion of the particle and the capillary 
waves leads to a rich behaviour. Before solving the Stokes 
equations, we first remark that the shape of the interface 
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depends on the detailed history of the motion of the par- 
ticle as well as on the shape at some earlier times. We 
are then naturally lead to perform a Fourier mode anal- 
ysis in time, the Fourier transform f(u>) of an arbitrary 
function f(t) being defined as 



dt exp[— icut]f(t) 



(5) 



Besides, one can note that the problem is translationally 
invariant along the direction parallel to the surface. It 
is thus helpful to use the two-dimensional Fourier repre- 
sentation introduced above in Eq. ©. It also appears 
judicious for this study to decompose the vector fields 
into their longitudinal, transverse and normal compo- 
nents [2^]. This defines a new set of orthogonal unit 
vectors (q, t, n), where q is the unit vector parallel to q, 
n the unit vector in the z-direction, and t the in-plane 
vector perpendicular to q and n. These vectors are con- 
nected to the cartesian basis (e x ,e y ,e z ) through 



e x + 



■2_ %_ _<h 



n 



(6) 



The velocity and the force are written v = viq+vtt + v z n 
and F = Fiq + Fit + F z n, respectively. Inserting these 
representations into the Stokes equations CQ) - 10 finally 
leads to a system of ordinary differential equations for 
the Fourier-transformed quantities 

W^ + laTJ ~iqp +F t 6(z - z ) = (7) 



-rj a q v t + r\ a 



dz 2 
d 2 v t 
dz 2 



+F t S(z - zo) = (8) 



-Vaq 2 v z 



1c 



dz 



, +F z 6(z~z Q ) = Q (9) 

dz 



with the divergenceless condition 

dv z _ 
dz 



iqvi 



. 



(10) 



Although this framework is not as "transparent" as the 
usual image method, its advantages are twofold. On the 
one hand, it is particularly well suited to accommodate 
with the description of the interface energy in Fourier 
space, since it thoroughly accounts for the symmetries of 
the problem. On the other hand, the transverse compo- 
nent of the velocity is decoupled from the longitudinal 
and normal directions. Moreover, relation (|10|l provides 
a usefull link between vi and v z , so that it is not difficult 
to get a single, fourth order differential equation for the 
normal component only 



d 4 v z 



,d 2 v z 



-q v z 



d(z-Zo)-\ 6 (z — zo) ■ 

(11) 



dz 4 dz 2 
Here, S' is the derivative of the Dirac delta-distribution 



III. BOUNDARY CONDITIONS 

To describe the flow in the presence of an interface, we 
must consider the flow on each side separately, and then 
require proper matching conditions for the velocity and 
surface forces. As stated above, the hypothesis of smooth 
deformations around the planar configuration enables us 
to express the boundary conditions at the undisplaced 
interface z — |20(- Because our representation of the 
velocity in terms of longitudinal and transverse coordi- 
nates is not commonly used in the litterature, we find it 
worthwhile to give some details regarding the derivation 
of the boundary values. 



A. Continuity of the velocity 

First of all, we have to ensure that the velocity is 
continuous at the interface. Explicitly, this requirement 
reads 

5i(q,0+ w)=«,(q»0".w) . (I 2 ) 
«t(q,Q+w) =«t(q,0-,w) , (13) 
v z (q.,0+,u) =v z {q,0~,cj) . (14) 

Interestingly, the condition (|12|) for the longitudinal coor- 
dinate together with the incompressibility condition 11(J|) 
implies an additional boundary condition for the normal 
coordinate of the velocity, namely 



dv z 
dz 



0+ 



dv z 
dz 



(15) 



B. Balance of tangential forces 



Secondly, tangential stresses have to be balanced at 
the interface. In real space, the continuity condition 
for the normal-tangential components of the stress ten- 
sor reads a zx \ + = <r zx \ - and a zy \ + = a zy \ -, with 
(Tjk = —pSjk + Va{dvj /dxk + dvk/dxj) the stress tensor 
in cartesian coordinates. Switching to {q, z, lu} variables, 
both requirements reduce to 



flvj| 

dz 



0+ 



= m 



<3vi| 
dz 



where the two-dimensional vector V11 = (vi,v t ) is the par- 
allel velocity. Projecting this equation onto the trans- 
verse direction leads to the condition 



dv t 
dz 



0+ 



dv t 
dz 



(16) 



whereas projection onto the longitudinal coordinate gives 
another condition which still involves both vi and v z . 
In order to obtain a boundary condition for the normal 
component only, the incompressibility condition (|10|l is 
once more invoked. We finally get 



dz 2 



+ q 2 v z 



= V2 



dz 2 



+ q 2 v z 



(17) 
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Note that the balance of tangential stresses is also rele- 
vant with regard to the normal component of the veloc- 
ity. 



C. Discontinuity of normal stress 

The next condition that has to be enforced concerns 
the normal-normal stress difference that comes into play 
whenever the interface is bent. Indeed, a deformation 
of the interface gives rise to normal restoring forces, 
expressed as the functional derivative of the Hamilto- 
nian J3J. For small displacements, the forces are small 
and proportionnal to h. The normal stress condition 
reads, in real space, a zz \o- — <J ZZ \$+ = —8H,/Sh. In terms 
of the variables {q, z, uj}, we have 



oz 



o+ 



n dv z 
oz 



= -E q h(q,uj) , 



where we define the energy density E q — j(q 2 + A -2 ). 
It can be seen that the normal stress difference at the 
interface is balanced by interfacial tension and buoyancy 
forces (due to the density difference between the two flu- 
ids). This condition still involves both the normal com- 
ponent of the velocity as well as the pressure field. To get 
a relation in terms of v z only, we shall first use Eq. Q) to 
express the pressure difference (remember that zq > 0) 



iq(p(0+)-p(0-)) 
'd 2 vi 



dz 



q 2 vi 



o+ 



dz 2 



q 2 vi 



Substitute vi for v z thanks to relation (|10|l . we arrive at 
the condition on the third derivative of the velocity 



d 3 v z 



-M 2 

d 3 v z 



,dv 2 



iq 2 



o+ 
dz 



-q 2 E q h{^uj). (18) 



D. Immiscibility of the two fluids 

To make the calculations tractable, we suppose that 
the condition of immiscibility can be written at z = 0. 
This approximation is justified since the fact that it is 
in any rigour valid at z = h is an effect of higher order. 
Within this assumption, the time rate of change of the 
shape function is related to the normal velocity at the 
interface through 



u z (q, 0, w) — iujh(q,u>) 



(19) 



up to linear order in the deformation field. This clo- 
sure relation is especially relevant since, as shown in the 
following, it allows to work out the instantaneous de- 
formation of the interface in response to hydrodynamic 
stresses. 



IV. GREEN FUNCTION AND 
TRANSLATIONAL MOBILITY 

A. Motion of the interface 

We now have all the ingredients to solve the Stokes 
equations. Because our calculations are algebraically in- 
volved, we save the details for the appendices. An inter- 
esting results is that the local deformation of the inter- 
face is directly proportional to the amplitude of the point 
force applied at height zq 



h(q,u) = R(q,z ,w)F(w) 



(20) 



where the vector R is the response function obtained 
thanks to the closure relation 1)190. For a vertical force 
F = (0, 0,F Z ), we find in Appendix El 



R z (q,z ,uj) = — 



1 



Anq (oj q + iuo) 



— (1 + qz )e" 



(21) 



As expected, the relaxation dynamics of the profile is 
governed by the mean viscosity rj — (r/i + r/2 ) / 2 . The 
response of a deformation mode with wavevector q is 
characterized by its relaxation rate 



(22) 



Remark that differents wavevectors are not damped in 
the same way. The amplitude of the response function 
is always maximum for q = 0, h(0,oj) = F z /(Apg). It 
then vanishes with increasing q, all the more rapidly as 
the frequency lo or the distance zq are large. The result 
Eq. (|21|l can be interpreted as follows. The real part of 
R z , which is in phase with the strain, is the analogous of 
a storage modulus for a viscoelastic medium [26^. This 
contribution corresponds to the elastic energy stored in 
the deformation of the interface. On the other hand, the 
imaginary part of R z plays the role of a loss modulus 
and describes the viscous dissipation associated with the 
relaxation of individual deformation modes. 

Coming back to the motion of the interface in real 
space, one can evaluate the inverse Fourier transform of 
the response function, though the calculations will not be 
performed here. A deformation may also be obtained as 
a result of a point force applied parallel to the interface. 
We find 



1 



Rl(q,z Q ,w) = - 

47]q [UJ — lUJq 

Rt(q,z ,uj) = , 



-qz e 



(23) 
(24) 



for the longitudinal and transverse coordinate, respec- 
tively. Note that the shape of the interface is not affected 
by the transverse component of the force. 
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B. Oseen tensor 



C. Translational mobility tensor 



The components of the Oseen tensor are then obtained 
by identification with the definition 



(25) 



where i,j G {l,t,z}. For symmetry reasons, the 
Green function satisfies the general relation Gij{r, r', t) — 
Qji(r',r,i) for z and z' on the same side of the 
interface This is re-expressed in terms of 

our particular choice of variables as £?y(q, z, z', ui) — 
Qji(—q,z',z,—w)*, property that can be checked along 
the calculations. As shown in the appendices, the Green 
function can always be written as 



Q(q,z,z ,u) 



^(q, Z -z )+Ag(q, Z ,z , W ) 
g^(q,z-z Q ) + Ag ( - 1 \q,z,z ) 
j(2) 



A£^(q,z,z ) 



(26) 



see Appendix ^ for the exact expression of Q t j, Ap- 
pendix[B]for the components Gzj , and Appendix^lfor the 
components Qij. The first term, which depends only the 
relative distance (z — zq), would reduce to the usual free- 
space Green function if the viscosities were equal. The 
second term, At/' 1 ', is the correction for an undistorted 
interface. Both contributions have already been obtained 
in previous work, though not in this particular choice of 
coordinates The original part of this study is the 

derivation the contribution coming from the deformation 
of the interface. The prefactor (w — iojq)^ 1 in Eq. 12(j|) 
is a clear signature of hydrodynamic scattering effects on 
the soft surface. Indeed, the fluid flow resulting from a 
displacement of the particle exerts stresses that deform 
the interface. Relaxing back to its equilibrium position, 
the interface creates a backflow that in turn perturbs the 
motion of the particle, and so forth. The infinite sum 
(u) — iujq)^ 1 = i/iOq Y^=q ( — *) n ( w / w <?)" i s the expression 
of this infinite series of "reflexions" of the original point 
force on the interface. This argument is confirmed by 
the fact that the correction vanishes for 7 — > 00. One 
recovers the results for the flat liquid-liquid interface in 
the high surface tension limit, as expected. 

Finally, once all the components are known in the 
(q, t, n) basis, it is not difficult to express the Oseen ten- 
sor in cartesian coordinates. In particular, the diagonal 
components are given by 



Gxx — -^Gu 



■>vv 



-^yu 



(27) 
(28) 



Similar relations can be deduced for off-diagonal terms, 
though they will not be required in the following. 



From the matrix elements of the Oseen tensor, we can 
obtain the mobility matrix for a sphere. To this aim, we 
still have to enforce the no-slip boundary condition for 
the fluid flow on the surface of the particle. In the follow- 
ing, we assume that the particle is a sphere of radius a. If 
we note U(ro) and ft respectively the translational and 
rotational velocity of the sphere, ro being the position of 
its center-of-mass, then the fluid velocity satisfies 



v(r + a) = U(r ) + ft x a 



(29) 



for any vector a scanning the surface of the bead. In- 
tegrating this relation over the particle, one obtains a 
linear relation between the total friction force F# ex- 
erted by the liquid and velocity of the particle |22| ■ This 
relation defines the (frequency-dependent) mobility ten- 
sor through U = —JlFh- It can be written as the sum 
of two terms, Jj, k i{zo,^) = Vofal + Apkl(zo,w), with 
fiQ = (Gnrjia)^ 1 the bulk value for a particle in fluid I 
but infinitely far from the interface. The correction A/Ijy 
is then expanded in powers of a/zo. In the limit of small 
particles a -C zq, the correction to the mobility tensor is 
given, at leading order, by 



A/Zfei (z q ,uj) 



(2tt): 



-A£ fci (q, z , z ,oj) . (30) 



As a matter of fact, it can be shown that all cross- 
contributions vanish, so that the correction to the mobil- 
ity tensor is also diagonal with elements Ap xx = Ajlyy = 
Ajuii and AJ1 ZZ = Ajuj_. 



Perpendicular mobility 



From the result Eq. (|B5|I for the normal-normal com- 
ponent of the Oseen tensor, we find 



A/ij_ (z ,u>) 



1 / 27?l+37? 2 

i67Tf7iz v m+m 
5 



£0 

lr. 



(31) 



In this expression, r = Arjl c /^j corresponds to the longest 
time required for elastic stuctures in the fluid — in 
our case, the interface — to relax. For typical values 
rj = I0~ 2 Pa.s and Ap = 10 2 kg.m~ 3 , it ranges from r s» 
10~ 3 s for usual interfaces with 7 = 100 mN.m" 1 up to 
twIs for ultra-soft interfaces with 7 = 0.1 /dM.m -1 |24| . 
The frequency-dependent contribution F arises from sur- 
face deformations and is therefore governed by the mean 
viscosity rj. It is given by 



F(s,k) 



dx- 



iksx 



r (l + kx) 2 exp[-2fcx] 



F'(s,k)+iF"(s,k) 



(32) 



6 



1 

0.8 

LL 

ffl 0.6 

CL- 
OT 

S 0.4 
0.2 




/ / 

/ / 
■■' / 
/ / 

y y — k=io 2 



k=0.1 

k=1 

k=10 



10° 1 o 2 

Reduced Frequency 



10 4 




Reduced Frequency 



FIG. 1: (Color online) Real part of F as a function of the 
reduced frequency s — lot, for different values of the reduced 
distance k = Zo/l c . The position of the cross-over between 
the two regimes lot <C 1 and lot ^> 1 is quite sensitive to the 
distance to the interface. 



This integral actually corresponds to the sum over all de- 
formation modes of the interface |2jj. For u> = 0, one has 
F(0, Zq/I c ) = and the authoritative reader will recog- 
nize on the right hand side of Eq. i|31|) the correction to 
the mobility of a sphere near a flat, liquid-liquid inter- 
face One even recovers the result of Lorentz for a 
hard wall by taking the limit r/2 — > 00. For finite values 
of lo, the additional term is actually a complex number. 
Its real part F' represents the contribution to the viscous 
dissipation that comes from interface deformations. As 
shown in Fig.^ F' is positive for any value of the param- 
eters, so that the real part of the mobility increases when 
the constraint on the shape of the interface is released. 
Viscous dissipation is therefore always lower for a soft 
interface, which can bend under hydrodynamic forces, 
than for a rigid interface. Another outcome is that the 
mobility of the particle also exhibits an imaginary part 
F", which corresponds to the storage of elastic energy in 
the deformation of the interface. As shown on Fig. [21 F" 
is non-zero only for intermediate values of the frequency 
lot ~ 1. The latter contribution vanishes when lo — > 00 
and one gets in this limit 



A/j,± (z ,lo -> 00) 



m - m 



i67r?7izo \m + m 



(33) 



Lastly we remark from Fig. ^ an d 121 that both F' and F" 
are of 0(1) for a wide range of reduced distances Zq/I c . 
But because of the prefactor z~ l in Eq. (|3*T)l . the coupling 
between the motion of the particle and the shape of the 
interface vanishes when the particle is far away from the 
surface, as one might expect. 

To better understand those results, let us focus on the 
relaxation time r. From the definition r = Arjlc/j <x 
it can be noticed that the limit lot <C 1 actually 
coincides with the limit 7 — > 00. At very low frequen- 
cies, the interface thus appears infinitely rigid and one 
therefore recovers the mobility of a particle near a flat 
interface. This also explains why the imaginary part of 



FIG. 2: (Color online) Imaginary part of F as a function 
of the reduced frequency s — lot, for different values of the 
reduced distance k — zq/I c . The elastic coupling is maximum 
around the value lot ~ 1. 



the correction vanishes when lo — 0. A similar reasoning 
applies to the other limit lot 3> 1, that would correspond 
to an interface with vanishing surface tension. Since no 
elastic energy can be stored anymore, F" has to vansih 
when lo — > 00. Interestingly, the sign of the real part 
of Apt^ may change depending on whether lot < 1 or 
lot 3> 1. Indeed, it is always negative at low frequencies, 
whereas it may be positive at high frequencies provided 
that 771 > r)2- This susprising property, peculiar to soft 
interfaces, may strongly influence the statistical proper- 
ties of Brownian particles since surface deformation may 
enhance diffusion - with regards to the bulk value - at 
short times. 



2. Parallel mobility 

Similar conclusions can be drawn for the mobility par- 
allel to the surface. We find that 



1 



2r/i - 3i] 2 
32irniz \ rji + ?/ 2 
1 



G[lot^ 
64^77^0 V 



(34) 



where the frequency-dependent contribution is given by 

iksx 



G(s,k) 



dx 



1 + isx + x 2 



fcV exp[-2fcr] . (35) 



In particular, one recovers the mobility coefficient for 
sphere near a rigid interface in the asymptotic limit 

LOT < 1 



A/xii {z ,lo = 0) = 



2?7i - 3772 



327r?7izo v m + m 

whereas one obtains in the other limit lot 3> 1 

A/i|| (z ,uj -> 00) = — ■ 

327r?7iZo V 771 +7? 2 



(36) 



(37) 
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V. DISCUSSION 

To summarize, we have calculated the mobility tensor 
of a spherical particle moving close to a fluid-fluid inter- 
face. Several lengths are inherent in the system, namely 
the radius a of the particule, the distance from the wall 
zq, and the capillary length l c . The results presented 
in this work concern the response to a point force and 
are valid for particles far from the interface a zq. Be- 
cause a soft interface can deform and store elastic energy, 
the mobility tensor decomposes into a real and an imag- 
inary part. In steady-state limit lot -C 1, deformations 
are irrelevant and one recovers the classical result for a 
flat, fluid-fluid interface. On the other hand, the short- 
time limit lot 3> 1 presents the intriguing feature that 
the perpendicular mobility can be higher than the bulk 
mobility if r\\ > 772. Yet this result does not break any 
fundamental law since it arises from the fact that the 
particle "feels" the other side of the interface, which has 
a lower shear viscosity. Finally coming back to time vari- 
able, the friction force experienced by the particle will be 
expressed as a convolution product and is therefore non- 
local in time. Solvent backflow and delay of the response 
of the elastic interface then induce memory effects in the 
motion of the particle. 

The framework developed in this study may be 
adapted to various problems near soft interfaces. For 
instance, one might investigate surface-mediated contri- 
butions to the coupled diffusion of two particles. One can 
also consider more complex surfaces, such as surfactant- 
covered interfaces or fluid membranes. Predictions re- 
garding the rotational mobility might be relevant for ex- 
periments as well, especially in the case of anisotropic 
particles. Remark that translational and rotational mo- 
tions are not coupled for a sphere in the linearized theory. 
This might not be true anymore for large deformations, 
where non-linear effects come into play |29| . 

Another point that might be included in the theory is 
the effect of fluid inertia. This contribution has been ne- 
glected so far, though it becomes relevant at frequencies 
higher than lo c = r//(pa 2 ). For typical values r\ = 1CP 3 
Pa.s, p = 10 3 kg.m -3 and a — 1 fim, we obtain io c « 10 6 
rad.s -1 . Here however, we consider time scales compa- 
rable to the relaxation time of the interface. This cor- 
responds to frequencies in the kH range, so that our ap- 
proximation is fully justified. At this point, it should be 
mentioned that a study similar to ours, including fluid 
inertia, has been published during the completion of this 
work 30]. In the steady limit, the author obtains the 
result for a rigid wall with stick boundary conditions. 
This however cannot be correct since one expects to find 
in this limit the mobility of a sphere near a fluid-fluid 
interface. The results derived in Ref. [3(| are therefore 
questionable, but a closer inspection would be required 
to identify the origin of the discrepancy. 

Finally, let us comment on some possible comparisons 
with experiments. Recently, de Villeneuve et al. have 
considered the sedimentation of PMMA spheres towards 



an interface with ultra- low tension 7 w 0.1 fiN/m [31j . 
In this regime, long-range hydrodynamic interactions are 
dominant and lubrication theory does not apply. The 
authors clearly observe strong deformations of the inter- 
face, of the order several micrometers for spheres with 
radius a = 15 /im |3lj . Moreover, they measure sedimen- 
tation velocities that do not follow the theoretical curves 
for an undistorted interface, the particles falling faster 
towards the soft interface. The interpretation of those 
results might be quite straightforward in the light of the 
present analysis, even though the non- linear equations 
of motion might be challenging to solve. Work on this 
question is currently under progress. 
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APPENDIX A: TRANSVERSE COMPONENT OF 
THE VELOCITY 

We begin with Eq. JSJ for the transverse component, 
which is easier to solve since it does not couple with the 
longitudinal and vertical coordinates of the velocity. This 
equation can be rewritten as 



d 2 v t 2 _ 
az z 



mi 



;S(z - Zq) 



(Al) 



With the condition that the fluid is at rest at infinity, the 
solution is 



vt(q,z,Lj) 



Ae~ qz 

Be qz + Ce~ qz 
De qz 



for < Zq < z , 
for < Z < Zq , 
for Z < < Zq . 



We then need to specify the boundary conditions in order 
to determine the four integration constants. The continu- 
ity of the velocity and the balance of tangential stresses 
at height z = give the conditions l|13|l and (|lfc>|) . We 
get another couple of conditions by invoking the stan- 
dard continuity conditions for the Green function at the 
singularity z = zq. Explicitcly, these requirements read 



dv t 
dz 



dv t 

dz 



V1Q 2 



(A2) 
(A3) 



Enforcing the boundary conditions (|13fl . (|16fl . <|A2I) 

and l|A3(l . we find the following expression for z > 



v t (q,z,Lo) 



F t 
2m q 



-g\z 



1 - A 
1 + A 



3 -q(z+z ) 



(A4) 
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and for z < 



F 



2rj 2 q \ l + \ 



,-g|*-*o| 



(A5) 



The transverse components of the Green function are 
then obtained by comparison of equation (|A4|) (for z > 
0) or (|A5|) (for z < 0) with the definition of the Oseen 
tensor v t = g tl F t + Q tt F t + G tz F z . Obviously, we get 
Gti = Qtz = 0, the only non-zero component being Gtt- 
Note that the transverse component of the velocity is not 
affected by the shape of the interface. 



APPENDIX B: NORMAL COMPONENT OF THE 
VELOCITY 

1. Differential equation and general solution 

The solution of the fourth-order differential equa- 
tion satisfied by v z is 

v z (q,z,u})=(A + Bz)e~ qz for z > z 

v z (q, z, w) = (C + Dz) e qz + (E + Fz) e~ qz for z < z 
v z (q, z,uj) = {G + Hz)e qz for z < 



For the sake of simplicity, we shall focus separately on the 
situations where (F; = 0,F Z ^ 0) and (F; ^ 0,F Z = 0). 
According to the superposition principle, each solution 
leads by identification to the components of the Green 
tensor Q zz and G z i , respectively. Obviously, the normal- 
transverse component is identically zero, G z t = 0. 



2. Normal-normal component 

We first consider the case where Fi = and F z ^ 0. 
In addition to the boundary conditions l|14|). I|15|). I|17|) 
and l|18|) . we have to enforce the usual conditions for 
the Green function at the singularity position z = zq, 
namely the velocity, its first and its second derivative are 
continuous at z = z a . The only discontinuity comes from 
the third derivative 



a 3 i 



dz 3 



d 3 v z 



dz 



(Bl) 



The algebra involved to evaluate the height integration 
constants is rather lenghty but presents no difficutly. We 
simply give the resulting velocity field 



v z (q,z,uj) 



F z 
\r) X q 



(1 + q\z - z |)e-^-^ - (j^jj i 1 + q{z + Z °) + 2q2zZ °) e- q{z+Zo) 



r 



u q h(<l,w)(l + qz)e qz , 
(B2) 



for z > 0, and 



F~ 



w g /i(q,w)(l - qz)e qz 



(B3) 



the velocity i|B2|) or (|B3|) at height z = and comparing 
with l|19(l then leads to 



/i(q,w) 



— L^(l + gz ) e -9*°ik. . (B4) 
ojq + 2lj Ar/q 



for z < 0. The velocity field still depends on the defor- Bringing Eq. (|B2|) and (|B3|I together with (|B4(1 . we fi- 

mation of the interface, which itself is a function of the nally obtain the normal-normal component of the Green 

velocity through the closure relation l|19|) . Evaluating function 

I 



G zz {q, z,z ,lu) 



for z > 0, and 



Szz(q, Z,Z ,Uj) 



1 



(1 + q\z - Zo\)e- q \*-*°\ -(l + q(z + Z ) + Y^f} e- q 



-q( Z + Zo ) 



+ 



1 w 



4rjq u) — iujq 



— (1 + qz)(l + qz )e- q ( z+zo '> 



1 



Ar] 2 q V 1 + A 



q 2 zz e q{z - za) + 



1 



4f]q u> — iu q 



(l-gz)(l + gz )e«( 2 - zo ) 



(B5) 



(B6) 



for z < 0. 
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3. Normal-longitudinal component 

In order to get the component Q n i of the Oseen tensor, 
we perform the same analysis expect that we now keep 
Fi 7^ 0, whereas we set F z = 0. This time, the disconti- 
nuity imposed by 8' in Eq. Ijllll has an incidence on the 
second derivative of the velocity at z — z 



dz 2 



d 2 v z 



(B7) 



the velocity, its first and its third derivative being con- 
tinuous. The algebra being quite similar to that of the 



previous section, we shall skip the details. Once again, 
the velocity field depends on the deformation of the inter- 
face. Intersetingly, a point force exerted parallel to the 
surface is responsible for a normal displacement of the 
fluid-fluid interface. Evaluating the velocity at height 
z = leads to the result 



/i(q, lo) 



to — lUJ a 



-qz e 



-qzo 



T=r ■ ( B8 ) 
Arjq 



Bringing everything together, we find the normal- 
longitudinal component of the Green function 



J 



£ 2 ;(q, z, z ,uj) 



q (z - z) e-^-^ + ( \-±qz - qz - ) ( ; 1 " 1 



1 + A 



1 + A 



Arjq oj — iuj q 



(1 + qz)qz e^ z+z °'> , 



for z > 0, and 



Gzl{<l,Z,Z ,L}) = -■ 



Ai] 2 q \l + X 



qz 



(1 - qz Q )e^ z - z ^ + 



Arjq oj — iuj q 



(1 - qz)qz e q{z - zo) 



for z < 0. 



(B9) 



(BIO) 



APPENDIX C: LONGITUDINAL COMPONENT pressibility condition JHJ}, vi is related to v z thanks to 
OF THE VELOCITY vi — (i/q)dv z /dz. From the definition of the Oseen ten- 

sor vi = Q u Fi + Q iz F z (since, of course, Q u = 0), it is 
To obtain the longitudinal component of the veloc- straightforward to get 
ity, there is actually no need to solve the correspond- 
ing differential equation Q. Indeed, from the incom- 



Giii^z^n,^) 



(l-q\z~z \)e-^ z - z °\ 



1 



1 - A 
1 + A 



l-q(z + z ) + ^^)r ".--..,> 



Arjq oj — iui q 



q zzoe 



-q(z+z ) 



for z > 0, and 



for z < 0. 



Gu(q, z,z ,uj 



Ar/q oj — iui q 



q zzoe 



q(z—zo) 



V 



(CI) 



(C2) 



Regarding the longitudinal-normal component, Gi z , no duced from Q z i using the symmetry relation of the Green 
additional algebra is required since it can directly be de- function. 
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